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We present the prospects for studying quantum correlated charm de- 
cays at the ?/>(3770) using 0.5-1.0 atT 1 of data at SuperB. The impact of 
studying such double tagged decays upon measurements in other charm 
environments will be discussed. 
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1 Introduction 



SuperB is a next-generation high-luminosity asymmetric-energy e + e~ collider that 
aims to collect 50 to 100 times more data than the B-Factories of the last decade, 
BABAR and Belle. The center-of-mass energy for the majority of the program will be 
at or near the T(4S) resonance with the designed peak luminosity of 10 36 s _1 cm -2 . 
The goal is to collect 75 ab _1 over five years. Additional runs are also planned at 
DD threshold ^(3770) to collect 0.5-1.0 atT 1 over a few months. SuperB's physics 
programs include, but not limited to, heavy-flavor B Uj d, s , D, and r physics. SuperB 
will be able to search for new physics at energy scale up to 10-100 TeV through 
rare/forbidden decay searches, CP violation, and precision CKM matrix measure- 
ments pp. 

With 75 ab -1 of data near T(AS), O(l0 n ) charm mesons will be created. They 
come from continuum production e~e~ — > cc, as well as B decays. Many charm 
analyses identify a D meson through D* + — > D°tc + process. Consequently the recon- 
struction efficiencies are relatively low. 

With 0.5 ab" 1 of data at ^(3770), one can expect approximately 1.8 x 10 9 D°D° 
and 1.5 x 10 9 D + D~ events. This amount is more than an order of magnitude 
larger than the current charm factory BESIII [4J will collect. As a J PC = 1 state, 
^(3770) — > DD is in a quantum entangled, anti-symmetric state. If one D decays to 
state a at time t\ and the other to /3 at £2, the decay amplitude M. is 

M = ^[{a\H\EP{t x )){p\H\W{t 2 )) - (P\H\D°(t 2 ))(a\H\D°(ti))]. (1) 

A neutral meson mixing system can be described by a 2 x 2 effective Hamiltonian 
with non-vanishing off-diagonal terms 

4( = p) 



The eigenstates ^1,2) = p\D°) ±q\D°) satisfy 



The eigenvalues are 

Ai, 2 = m 1>2 - l -T 1>2 =(^M- l -Y^± q - (m 12 - l -T 12 ^j . (4) 

Here we have assumed CPT conservation. The time evolution of ^(3770) — > DD — >■ 
a(ti)(3(t 2 ) system can then be expressed as 

dT/dt oc + | 2 + |a_ | 2 ) cosher At) + (\a + \ 2 - |a_ | 2 ) cos(xrAt) 

-2TZe(a* a_) sinh(yrAt) + 2Xm(a*_a_) sinh(a;rAt), (5) 
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where At = t 2 -h, a + = A a A -A a A fi , a_ = -jA a A^ + ^A a A , M = (M n + M 22 )/2 
r = (Tn + r 22 )A x = (mi - m 2 )/r, and y = (T 1 - T 2 )/(2T); A X (A X ) is the decay 
amplitude of D (D) to X. 



2 Charm mixing measurements 

The mixing in neutral D system is expected to be very small. The short- distance 
|AF| = 2 comes from box diagrams. The diagrams with b quark in the loop is 
CKM-suppressed (with V u b in the vertex), and the ones with s and d quarks are 
GIM-suppressed. The long-distance contributions come from diagrams that connect 
D° and D° through on-shell states (e.g., KK). These long-distance contributions are 
expected to be 0(lCr 3 ) but the theoretical calculation is difficult [3]. CP violation 
induced by mixing is therefore expected very small too. Observations of large mixing 
and/or CP violation are considered clear signs of new physics beyond the standard 
model. 

Charm mixing has been firmly established at 5-factories [2] using continuum 
events in the data taken near the Y(4S) resonance. Both x and y terms are approxi- 
mately 0.5%. These analyses use the charge of the soft pion from D* + — > D°tt + (or 
its charge conjugate process) to identify the initial flavor of the D meson, and recon- 
struct a final state that is accessible by both D° and D°. The decay time distribution 
for the D meson tagged by tt + (N(t)) and n~ (N(t)) are 



N(t) oc 
N(t) oc 



■>■- ■ >r , |2rpjl 2 



i + —£*-\\ f \*(rty + \x f \( y cos(5/ + <p f ) -xsin(s f + 4> f )(rt) 

2 2 

1 + X -Y-\\ f \'\Ttf + |A / r 1 ( 2/ cos(^07) - xsin(^07)(rt) 



(6) 
(7) 



where xTt, yTt <C 1, A/ = (qAf)/(pAf), and 4>f — i>f + 4> m , where 5f (i/jf) is the 
relative strong (weak) phase in decay, and m is the mixing phase aig(q/p). 

The cleanest mode used in this method is D — > K + it~ , which is Cabibbo-favored 
in D° decays but doubly-Cabibbo suppressed in D° decays. One does not measure 
x and y directly. Rather, the observables are rotated by the strong phase difference 

x' = x cos 5k-k + y sin Sk-k (8) 
y' = ycos5 Kn - xsmS K7T . (9) 

Independent measurements of strong phase difference are needed. 

Strong phase differences can be measured in ?/>(3770) — > DD decays with a 
"double-tag" technique. Due to the quantum-entangled nature of the system, when 
one D decays to a CP final state, the other D is projected to the orthogonal state, 
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which is a linear combination of D° and D°, and its decay branching fraction is sensi- 
tive to the relative strong phase of D° — > f and D° — > f. For example, for / = K~tt + , 
the effective branching fraction of the double-tag event is [5], 

J r s ± ,K-n+ - B s± B K - n +(l±2r cos 5 Kn + R ws + y), (10) 

where B$ ± and Bk--k+ are the branching fractions of D° decaying to CP± and K~tt + 
final states, respectively, (K + n~\D ) / {K + n~\D ) = re~ lSK7T , and Rm is the wrong- 
sign total decay rate ratio, R M = T(D° K--k + )/Y(D° K~tt + ) = r 2 + ry' + 
(x 2 + y 2 )/2. CLEO-c [5] has demonstrated this technique with 281 pb~ x of data and 
obtained 5 Kw = (22±^ii)° or [-7°, +61°] interval at 95% confidence level. 

Another powerful method of measuring D°-D° mixing is using a time-dependent 
Dalitz-plot analysis with three-body decays. With this method, one can avoid strong 
phase ambiguity and resolve x and y by exploiting strong phase variation and in- 
terferences of resonances on the Dalitz plot. The most power mode of this kind is 
D° — > K®n + ir~ . The time-dependent decay amplitude of a state created as D° or D° 
at t = can be expressed as [BJ, 

M(s 12 ,s 13 ,t) = A D (s 12 ,s„ f^ + + g -A D (s l2 ,s 13 f' {t) ~ e ^\ (11) 
M(s 12 , Sl „t) = A D ( Sl2 , Sa f^ + »(«) + lA D ( Sl2 , Sl3 f>V - e *\ (12) 

where Ad (Ad) is the decay amplitude of D° (D°) as a function of invariant mass 
squared s 12 = m 2 : = (p K o + p w -) 2 , s 13 = m\ = (p K o + p w +) 2 , and e 1>2 (t) = 
exp[— i(rrii t 2 — ir 12 /2)t]. Using this method, Belle [7] and BABAR [5] measured 
x = (0.80±0.29l°$i°;^)%, y = (0.33 ± 0.24±g:^±g;g|) % , andx = (0.16±0.23±0.12± 
0.8)%, y = (0.57 ±0.20 ±0.13 ±0.07)%, respectively, where the first uncertainties are 
statistical, the second are systematic, and the third are Dalitz plot model uncertainty. 

With 75 ab -1 of at T(AS) at SuperB, the statistical uncertainty can be reduced by 
a factor of 10. Since major systematic uncertainties are in fact statistical in nature, 
estimated from data control samples and simulated events, they will also be improved 
with more data. However, the Dalitz plot model uncertainty may not improve much 
without other input; it will become the dominant uncertainty at SuperB pp. 

To avoid Dalitz plot model dependence, Giri et al [9] proposed a method, origi- 
nally for measuring the CKM angle 7 in B + — > D[K®7t + tt~]K + decays using time- 
dependent Dalitz plot analysis. In this method, the Dalitz plot phase space is divided 
into N pairs of bins; two bins in each pair is mirror-symmetric about the line su = S13 
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the Dalitz plane. One then can define 

Cj = / _ dpAi2,i 3 ^13,12 COs(5i2,i3 ~ #13,12), (13) 

«i = J dp Ai2,i3A 13tl2 sm(5 12 ,i3- 813,12), (14) 

G?M?2,13, (15) 

where #ij,ifc = 5(sij,sik), and Ai^i* is the magnitude of the D decay amplitude 
A D (sij, Sik) = Aij t \k exp(z5y )1 fc). The integral is over the phase space of the bin i. 
Here we have used the fact that Ad(si 2 ,si 3 ) = Ao(si3, s 12 ). The Cj and contain 
unknown strong phase difference £12,13 — #13,12, and thus unknown, but T, can be 
measured with flavor tagged D° decays. For mirror bins, i and i, Ci = cj and = — sj. 
With charm mixing, the number of events in bin i at time t is [10] 



7?(*) ex e- r '[T, + VWiiciy + s lX )Tt + 0((x 2 + y 2 )(Tt) 2 )]. (16) 

One can fit all bins simultaneously to extract mixing parameters (x,y) if (sj,q) are 
known. 

Again, using entangled ^(3770) — > DD, one can measure and q. If one D 
decays into a CP eigenstate, the other D is in an orthogonal state. We denote these 
two states as D± = (D° ± -D )/ a/(2). The amplitude and partial decay width of the 
second D can be written as [9] 

A(D° ± ir s °(pi)7r-(p2)7r + (p 3 )) = ^=(^( Sl2) s 13 ) ± A D ( Sl3 , Sl2 )), 

dT(D° ->K s °(pi)7r-(p 2 )vr + (p 3 )) = 

X 

2 (^12,13 + ^13,12) ± ^12,13^13,12 COs(5i2,13 ~ 5 13 , 12 )dp, (17) 

where p$ in parentheses are the momentum of the corresponding particle. We can 
then measure Cj using 



JdT{D% ^ s ( Pl )7r-(p 2 )7r+(p 3 )) - JdT(D°_ ^(piKfeMPa)) 



(18) 

If we can bin the Dalitz plot so that q and Sj are nearly constant in each bin, (c^ Sj) 
can be determined with high precision 

Ci = Cj = cos(Sj — 6j)Apj = ^TfT-^ cos(#j — <5j), (19) 



* = J]^7)sin(5 J -5 J ) = ^± v / T J T J -c J 2 . (20) 
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Table 1: D°-D° mixing parameters (x,y) and strong phases obtained from x 2 fits to 
observables obtained either from BABAR or from their projections to SuperB. Fit a) 
is for 482fb _1 from BABAR alone and this is scaled up in b) to 75 ab -1 at T(4S') for 
SuperB. Fit c) includes strong phase information projected to come from a BES III 
run at DD threshold, and d) is what would be possible from a 500 fb -1 DD threshold 
run at SuperB. The uncertainties due to statistical limitation alone are shown below 
each fit result. 



CLEO-c pU 12] measured and a for D K° s rr + 7r~ and D K S K + K~ with 
818 pb~ x of data on -^(3770) resonance. They also estimated the impact on the 
measurement of the CKM angle 7. They found their Sj and q are consistent with 
that calculated from the Dalitz plot model used in BABAR analysis, and the reduction 
of Dalitz plot model dependence is substantial. 

3 Projected precisions in SuperB era 

In SuperB's physics reach studies [TJ, the expected precisions in D°-D° mixing param- 
eters in various scenarios on the SuperB time scale are estimated. First the results 
from BABAR's 482 fb -1 are extrapolated to SuperB's target of 75 ab -1 near T(4S"), 
without any independent inputs of strong phase measurements. Then the improve- 
ment due to better precision in strong phase measurements using DD threshold data 
are estimated, first from the forthcoming BESIII runs and from SuperB plan (0.5 ab -1 
integrated luminosity) . 

The results, including current average values from BABAR, are summarize in Ta- 
ble [TJ and the corresponding confidence regions are shown in Fig. [IJ 
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Figure 1: The confidence regions of D°-D° mixing parameters (x,y) in various 
scenarios described in the text and in Table [H Shaded areas indicate the cover- 
age of measured observables lying within their 68.3% confidence region. Contours 
enclosing 68.3% (la), 95.45% (2a), 99.73% (3a), 99.994% (4a) and 1 - 5.7 x KT 7 
two-dimensional confidence regions from the \ 2 fit to these results are drawn as solid 
lines. 

4 Time-dependent CP asymmetry 

Using coherent ^(3770) — > D°D° decays, one can perform time- dependent CP asym- 
metry studies analogous to T(4S) — > B°B° in 5-factories. If a neutral D meson 
decays to a final state at t\ that can identify the sign of its c-quark, e.g., lepton 
charge in semileptonic decays, the other D meson must be in an orthogonal state, 
i.e., the opposite flavor to the first D. The time-dependent decay rate of the second 
D meson into a CP eigenstate can be derived from Eq. |5j 

_ T(At) - r(At) _ rAt (\X f \ 2 - 1) cos(aTAt) + 21mX f smjxTAt) 
{ 1 F(At)+r(At) (l + |A / | 2 )(l + e 2 ^ rA *) + 2{1- e*y™)Ke\ f ' 1 1 
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Figure 2: Distribution of time-dependent CP asymmetry for D° mesons with TZeXf = 
ImXf = 1/V2 and CP = —1. These distributions assume q/p — 1 and the solid 
(dashed) line corresponds to y = 0.01(0.00). 

where At = t 2 — ti, and A/ = (qAf)/(pAf). 

Measuring time-dependent CP asymmetry in D°-D° system is much more difficult 
than in B°-B° system. The reason is that charm mixing rate is very small; both x 
and y are (9(1%) for D°-D°, whereas x ~ 0(1) for B°-B°. This effect is illustrated in 
Fig. [2] [13], in which one can see that even with a large CP- violating phase (arg(A/) = 
7r/4) the CP asymmetry is only a few percent within \At\ < 10 ps (more than 20 
times the D° lifetime). In contrast, within the same At range, the CP asymmetry for 
B° meson exhibits 1.5 full sinusoidal oscillations already. 

At SuperB, the design beam spot is much smaller (a x ~ 8 /im, a y ~ 40 nm, 
a z ~ 200 /im) than that in BABAR. This makes fitting for the primary vertex possible 
(and meaningful) even if no charged tracks originating from the primary vertex. As 
illustrated in Fig. [3j the charm mesons from ^(3770) decay fly away from the primary 
vertex for (9(100 /^m), depending on the center-of-mass frame boost. One can perform 
a beam-spot constraint fit on the ■0(3770) — > DD system to simultaneously fit for 
both flight lengths (L\ and L2) and convert them to decay times. At near T(AS), one 
studies charm physics using continuum data; the soft pion from D* decays are used 
to identify the initial flavor of the charm meson. 

SuperB has conducted studies to evaluate the sensitivities to mixing parameters 
(x,y) and CP-violating parameters q/p using 0.5 ab -1 of DD threshold data alone 
and compared that with using 75 ab~ x of data near T(AS). The preliminary results 
that used several two-body charm decays with various combination of CP/flavor-tags 
shows that the uncertainties with ?/>(3770) data are about six times larger than those 
with T(4S*) data. On this topic along, 0.5ab _1 of DD threshold data is clearly not as 
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Figure 3: Illustrations of charm meson reconstructions with beam spot at SuperB for 
(left) ^(3770) DD events and for (right) continuum events near T(4S). 

competitive as T(4S) data. However, one should be reminded that the former only 
requires a few months of data taking, while the latter will take five years according 
to the plan. 



5 Summary 

The precision of charm mixing measurements will be limited by the uncertainties of 
strong phases and Dalitz plot model by the time SuperB collected its targeted data 
near T(4S). One can mitigate this situation by utilizing the quantum correlation 
of charm decays in ^(3770) — > DD with BESIII data. With a months-long run at 
DD threshold at SuperB, it is possible to improve the precision by another factor 
of two. With a boost of the center-of-mass frame, time-dependent CP asymmetry 
measurements can also be performed in ?/>(3770) — > DD data, but the precision is not 
as competitive as the much larger T(4S) data. 

Finally, not discussed in this paper but worth noting here, charm threshold data 
have advantages to T(4S) data in several areas due to the low background and the 
fact that the whole event can be fully reconstructed (double tag), in addition to 
the quantum correlation. These areas include rare decays (D° — > 77, /z/i(X), etc.), 
leptonic/semileptonic charm decays, form factor measurements, CPT violation, CP 
violation in D — > that probes chromomagnetic dipole operator [jj. ] . and others. 
It certainly adds to the breadth of SuperB physics programs. 
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